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We propose an alternative notion of time reversal in open quantum systems as represented by 
linear quantum operations, and a related generalization of classical entropy production in the envi¬ 
ronment. This functional is the ratio of the probability to observe a transition between two states 
under the forward and the time reversed dynamics, and leads, as in the classical case, to fluctua¬ 
tion relations as tautological identities. As in classical dynamics in contact with a heat bath, time 
reversal is not unique, and we discuss several possibilities. For any bistochastic map its dual map 
preserves the trace and describes a legitimate dynamics reversed in time, in that case the entropy 
production in the environment vanishes. For a generic stochastic map we construct a simple quan¬ 
tum operation which can be interpreted as a time reversal. For instance, the decaying channel, 
which sends the excited state into the ground state with a certain probability, can be reversed into 
the channel transforming the ground state into the excited state with the same probability. 

PACS numbers: 05.30.d, 03.65.Ta, 05.40.a, 05.70.Ln 


Introduction: The discovery of fluctuation relations [l]- 
@] has transformed classical out-of-equilibrium thermody¬ 
namics, giving rise to the new field of stochastic ther¬ 
modynamics [Jt 7|. Important results obtained in the 
last two decades include the use of Jarzynski equality 
to measure equilibrium free energy differences in large 
biomolecules from non-equilibrium measurements [§, ij], 
generalizations of fluctuation-dissipation theorems from 
the equilibrium to the non-equilibrium domain li, 
and a sharpening of Landauer principle on the minimal 
heat generated in computing jl2Hl4j . The central quan¬ 
tity in stochastic thermodynamics is the entropy produc¬ 
tion in the environment, a functional of the whole system 
history ii which can be defined in two ways. The 
first method is by Clausius’ relation 

SS env = (3SQ (1) 

where /3 = j s the inverse temperature and 5Q is the 
heat (111 and the second way is as the Radon-Nikodym 
derivative of a forward and a reversed path probability 

IS,„„ bath] = log Ag, (2) 

Here P F [path] is the probability of the forward path and 
P fi [path] is the probability of the time-reversed path 0 , 
and fluctuation relations follow from © as mathemati¬ 
cal “tautologies” [l4], S3- Physically, fluctuation rela¬ 
tions are, however, not tautologies, because the quanti¬ 
ties in © and © should be the same. For standard 
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Markov models of the system-bath interactions (master 
equations, diffusion equations), this is indeed the case, 
but for more general models the situation is less evident. 

The (possible) extension of fluctuation relations to the 
quantum domain has been the focus of intense investi¬ 
gations reviewed in 0i3. Except for the generaliza¬ 
tions of the Jarzynski equality and Crooks’ fluctuation 
theorem to closed quantum systems [2l[ the results ob¬ 
tained to date lack the generality and simplicity of fluc¬ 
tuation relations in classical systems, and typically hold 
for specific models such as e.g. when the quantum jump 
method [22l l2jl | or the Lindblad formalism [2j| can be 
applied. 

In this work we focus on the tautological aspect of 
quantum fluctuation relations, i.e. on the analogues of 
©, which have not, we believe, been sufficiently em¬ 
phasized in the literature. To do this we have to define 
a general notion of time reversal of open quantum sys¬ 
tems. As in the case of a classical system interacting 
with a heat bath this notion of time reversal, which we 
call the R operation, is not unique. In the sense intro¬ 
duced here an open quantum system can be time reversed 
in many ways [261 ], for each one one can define an entropy 
production functional analogous to © and obtain fluctu¬ 
ation relations. We first give a very general (permissive) 
definition of the R operation, and show that it always 
leads to fluctuation relations. We then turn to a possi¬ 
ble definition of R starting with the standard quantum 
mechanical time inversion of a combined system and the 
bath and continuing to intrinsic representations in terms 
of Kraus operators. We discuss several special cases such 
as unital maps, and give examples of time-reversals of 
1-qubit channels. 

Generalities and definition of R operation: A state of 
a quantum system is described by its density matrix p 
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which is an TV-dimensional positive Hermitian operator 
with unit trace. Any physical operation on a quantum 
system can be described by a completely positive linear 
map which preserves the trace and sends one den¬ 
sity matrix p into another state p' = <b(,o). We define a 
general time reversal R as an involution on the set 12 n 
of quantum operations i.e. a bijective transformation 
R : $ —y <h R which satisfies (<& R ) R = $. 

R operation and fluctuation relations: To arrive at 
fluctuation relations we consider the paradigmatic exam¬ 
ple of two measurements, one before the beginning of 
the process described by 4>, and another after (2lj| . We 
will denote by A and O two measured operators (with 
eigenstates | a) and |o) and eigenvalues E a and E 0 ). At 
the beginning of the process the system will hence be 
in the pure state |a)(a|, and the probability to observe 
o at the end will be (o|<I>(|a)(a|)|o). For the reversed 
chain of events, we first measure O to obtain o, then act 
on the system with the reversed quantum map and 
measure A to obtain a. This happens with probability 
(a|$ fl: (|o)(o|)|a). A straightforward generalization of © 
to the quantum domain is thus 


dS en v \o: o] 


( 0 |$(|q)( 0 |)| 0 ) 

(a|<f>- R (|o)(o|)|a) 


( 3 ) 


To derive fluctuation relations from m one proceeds 
by analogy with the classical case using ©. For 
Jarzynski equality one takes A as an initial Hamiltonian 
Hi and O as a final Hamiltonian Hf with thermody¬ 
namic equilibrium states pf and pj, respectively, and 
defines the work done on the system during the pro¬ 
cess as 5W[a,o] = (E a — E a ) + /3~ 1 SS env [a,o\. Then 
0 e-P E °{a\<& R (\o){o\)\a) = Zf and this identity can 
simply be rewritten as 


]T(a|pf |a)(o|$(|a}(a|)|o)e-^-“’°e-^^-^) = e~^ F . 

a,o 


Here A F denotes the difference between the free en¬ 
ergy of both thermal equilibrium states, and the left 
hand side can be interpreted as (e~ ljSW ). For Crooks’ 
relation one similarly defines P F (x ) = (5(6W — x)) eq 
and P R (x) = (6 (6W R — x)) R q where in the latter the 
time-reversed work is 8W R [o,a] = -W[a,o] and the 
average is performed over the final equilibrium state 
pj and the reversed quantum map This implies 

P R (—x) = e~P( x ~ AF )P F (x), which is the relation of 
Crooks 

Standard quantum mechanical time inversion and time 
reversal of environmental representations Consider first 
a closed system. If p' = $p = U pW where U = 
f~e~'a^ tf ~ ti ^ H is unitary development with Hamiltonian 
H then the final state is time inverted by an antiuni¬ 
tary operator 0 i.e. p' = 0p'0 _1 . Acting on this state 
with U = 0LU0 -1 , we obtain a state at the initial time 
p = Up'U' and time inversion is complete if p = 0 _1 p0, 
which is the case. The time reversal of a closed system, 


based on standard quantum mechanical time inversion, is 
then the operation <f> —» <f> where <f>p = UpUf Obviously 
this is an involution: performing it twice we get back the 
map <f>p. 

Quantum maps have (generally non-unique) environ¬ 
mental representations 

p' = $p = Tr B [U(p®a)Ul] (4) 

where the principal state p acts in the space Ha while 
the ancillary state cr acts in the space Hb describing the 
environment. Both subsystems are coupled by a unitary 
operation 17, and the image p' is obtained by tracing out 
the environment. Under standard quantum mechanical 
time inversion U transforms to U and cr to a which allows 
us to define an environmental time reversal of the open 
quantum system as 


$ R Vp = Tr B 


u (p® a) m 


( 5 ) 


The disadvantages of such a definition are obvious: it is 
contingent on the chosen environmental representation, 
and to arrive at an intrinsic notion we have to pursue an¬ 
other approach. In addition, also when an open quantum 
operation has a very natural environmental representa¬ 
tion, © is not the only natural notion of generalized time 
reversal. We will return to this point below. 

Kraus form of quantum maps and the dual map: In¬ 
stead of © we will now start from the Kraus form (28j 
of quantum operator 4>: 

k 

P 1 = $(p) = J2 A iP A l ( 6 ) 

1=1 

As a quantum operation <f> preserves the trace, Trp' = 
Trp, the Kraus operators Ai satisfy the identity resolu¬ 
tion AJ A, = 1. The trace preserving conditions in¬ 
duce N 2 constraints so the set 12 a? of quantum operations 
acting on TV dimensional states has TV 4 — TV 2 dimensions. 
Although the number k in © is arbitrary, for any map 
there exist the canonical Kraus form , for which all Kraus 
operators are orthogonal, 

Tr A i A]=d i S i , j , (7) 

so that the Kraus rank k < TV 2 . In a generic case the 
numbers di are different and this representation is unique 
up to the choice of the overall phases of the Kraus oper¬ 
ators - see e.g. [29| . 

In the case of a unitary evolution, p' = ^u{p) = UpUf 
one has k = 1 and the only Kraus operator is unitary, 
A\ = U. Unitary maps belong to broader class of unital 
maps , which preserve the identity, 42(1) = 1, and satisfy 
the dual condition ]T7=i A t A\ = 1. A map which is trace 
preserving and unital is called bistochastic, as it forms a 
quantum analogue of a bistochastic matrix, which acts 
in the set of TV-point probability vectors. 

For any quantum map $ in the form © one defines 
its dual map 4>° such that 4>°(p) = A\pAi. Writing 
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for short $ = {Ai,..., A/-} we have $° = {A \,..., A[.}. 
Note that if a map $ preserves the trace, the dual map 
<h° is unital. If $ is bistochastic, so is <f>°. It is also 
convenient to define the class of self dual maps which sat¬ 
isfy <f>° = <f>. A quantum operation for which all Kraus 
operators are hermitian is bistochastic and selfdual. So 
is a map for which all non-hermitian operators occur in 
pairs, e.g. = {A ± , .. .,A m ,A\, ..., A^}. 

Restricting our attention to bistochastic maps as a first 
generalization of the standard time inversion, we can de¬ 
fine a time reversal Q Rb := <f>°. The choice Rb implies 
that (o|$(|<z)(a|)|o) = \{o\Ai\o)\ 2 = (a|$«(| 0 )(o|)|a). 

Therefore the fraction in Eq. m is equal to unity 
and hence for any bistochastic map $ and any initial 
and final pure states the entropy production vanishes, 
SS env [a, o] = 0. Although bistochastic maps are not time 
invertible in the standard sense, they are therefore time 
reversible in the sense that there exists a (natural) time 
reversal of such maps for which the entropy production 
in the environment vanishes. As a consequence fluctua¬ 
tion relations hold for all such maps with the work taken 
equal to the internal energy change which reproduces the 
original quantum fluctuation relation of Kurchan for uni¬ 
tary maps 0 as well as the result of Rastegin for bis¬ 
tochastic ma ps |3C| . In this connection it was recently 
observed |3l], [HI that the Jarzynski relation, with work 
taken equal to internal energy change, in general does not 
hold for nonunital maps. Thus if the map is nonunital 
then its dual is not a well-defined quantum map, time 
reversal cannot be defined as a dual, and the ratio in Eq. 
d3]) may be different from unity. 

A more general definition of the time reversal operation 
for a non-unital a quantum operation $ = {Ai,..., A k } 
was proposed by Crooks (26j. It is based on the invari¬ 
ant state of the map, p * = $(p*). The reversed map 
is given by the following sequence of the Kraus opera¬ 
tors $ Rc = {Af,...,A£}, where A? = pl /2 A\p~ 1/2 . 
Note that in the case of unital maps one has pj 2 = 1/TV 
so that Af = A] and one arrives at the dual map, 
= <j)° However, this definition does not work, if 
the invariant state belongs to the boundary of the set of 
quantum states so that p* is not invertible. 

The isomorphism of Choi and Jamiolkowski fiol liidl. [if i| 
states that an operation $ can be uniquely described by 
a dynamical matrix, (Choi matrix) 

D $ = TV($® J)|V’ + )(V’ + |, 

where \ip + ) = Y^iLi I*) ® |«) denotes the maximally 
entangled state on an extended Hilbert space Ra ®Rb- 
The matrix D of order TV 2 is hermitian and positive by 
construction, and its eigenvalues determine the relative 
weights di , while the eigenvectors of length TV 2 , reshaped 
into square matrices of order TV and rescaled by \fdt yield 
the Kraus operators Ai in the canonical form 0. 

Two maps and $2 are called unitarily equivalent, 
written dq ~ $ 2 , if there exist two unitary matrices V) 
and V 2 such that $ 2 (p) = V 2 {^\ (VipV^)) so the map 


<f >2 can be written as a concatenation of $1 with two 
unitary operations, $2 = 0 $1 ° ^Vi . Observe that 

for any two unitarily equivalent maps the corresponding 
dynamical matrices, D\ and D 2 , share the same spectrum 
and are unitarily similar. 

Essential map and its time reversal: We now look for 
a possible choice of the involution R for a general, non¬ 
unital quantum map, for which <t>° </ Qn Consider a 
generic quantum operation $, for which the spectrum 
{di }\Li °f the corresponding dynamical matrix is 
non-degenerate. The trace of the dynamical matrix is 
fixed, TtD = 2_, i=1 di, so let us order the Kraus oper¬ 
ators forming the canonical form 0 according to their 
norms, di = ||Ai|| 2 > d 2 = \\A 2 \\ 2 d k = \\A k \\ 2 . 

The leading Kraus operator, with the largest norm, is 
represented by a possibly non-hermitian matrix A± of or¬ 
der TV, which can be brought to the diagonal form by the 
singular value decomposition, 

Ai = Hi EV$, where ViV? = V 2 V^ = 1. (8) 

Here E is a diagonal matrix with all non-negative entries. 
In the generic case the spectrum of E is non-degenerate, 
and this decomposition is unique up to the phases of 
the right and left eigenvectors of A\ which form unitary 
matrices V\ and V 2 . In the degenerate case this decom¬ 
position is not unique. For instance, if A\ = U is unitary, 
than E —t, and one can choose e.g. V\ = U and V 2 = 1. 

For any map $ we select in this way two unitary ma¬ 
trices Vi and V 2 , which allow us to define rotated Kraus 
operators Bi = V^AiV 2 and the essential map 

k 

Hp) = B iP B\ = V^A t V 2 pV^A\Vi. ( 9 ) 

i—i 

For any operation determined by a set of ordered Kraus 
operators, $ = {Ai, A 2 ,..., A k } the corresponding es¬ 
sential map reads thus $ = {E,V^ A 2 V 2 ,... ,V^ A k V 2 }, 
where A\ = V\ EV 2 ■ Observe that the map $ and its 
essential map are unitarily equivalent, 

$ ~ 4> = o $ o . (10) 

It is easy to see that for any unitary evolution >F u the 
corresponding essential map reduces to identity map, 
\k (7 = 1. Thus the essential map generically provides 
a unique description of non-unitary part of a discrete 
quantum evolution. 

Consider, for instance a one-qubit Pauli channel, 

4 

®p(p) = 

3 =1 

where pj is an ordered, p\ > p 2 > p% > P 4 , normalized 
probability vector, ]Oj=i Pi = 1, while hermitian Kraus 
operators Aj form an arbitrary sequence of three Pauli 
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matrices cr*, i = 1,2,3 and identity, ero = 12 - Then the 
corresponding essential map reads 

3 

$p(p) = POP + ^2 WiPVi > ( 12 ) 

1=1 

where the three coefficients qi are up to a permutation 
equal up to the three smallest coefficients P2,P3,Pa of the 
map <I>. 

In other words any Pauli channel can be represented 
by a point in the regular simplex A 3 C R 3 spanned by 
the identity and Pauli matrices (see Fig. la). The corre¬ 
sponding essential map belongs then to the asymmetric 
fourth part of the simplex with a corner representing the 
identity map 1 . 



FIG. 1: a) The set of one-qubit Pauli channels forms a prob¬ 
ability simplex A 3 spanned by identity and Pauli matrices, 
b) An asymmetric part of the simplex forms a fragment of 
the set of essential maps 4>j, which contains the identity map 
1 and the maximally depolarizing channel <f>*. Channel A xy 
forms the symmetric combination of o x and a y . 

The essential map and an intrinsic R: Having defined 
an essential map $ , which describes the non-unitary part 
of the evolution, $ = 'F^t o 4> o 'I'y 2 , we are in position 
to give an intrinsic definition of a time reversed map 

$ R ■- \& V2 o $ o $ t . (13) 

Looking at the reversed map in the Kraus representa¬ 
tion = {A\, YA 2 Y.... YAkY} with the unitary ma¬ 
trix Y = V 2 V 1 , we see that the leading Kraus operator 
A\ is indeed inverted into A\, while other operators are 
suitably rotated to keep the map <F i? trace preserving. 

Making use of the form ca we see that the composi¬ 
tion of a map with its reverse reads 

o $ = ^v 2 ofolo vj/^t. (14) 

and is unitarily similar to <3E> o $ = $ 2 . 

It is easy to check that for any map the following prop¬ 
erties hold true ($ fl ) = <l> = (4>) fl and (<b R ) R = $, 
so that the reverse operation R is an involution, as re¬ 
quested. If the map is unitary, $ = iky, such a composi¬ 
tion reduces to the identity map 

= ( 15 ) 


as any unitary evolution can be reversed. 

Note that for a selfdual map one has = <f>° = $ 
so that = <f> 2 . For instance, consider the selfdual 

maximally depolarizing channel, which sends any state 
into maximally mixed state, <F*(cr) = 1/A. For any two 
pure states one has |(f|^ > *(|j)(i|)H)| = 1 /A so that the 
fraction in Eq. m is equal to unity for any choice of 
initial and final states. 

For any stochastic map described by two Kraus op¬ 
erators, $ = {Ai,A 2 }, the reversed map reads = 

{A \, — AiA\}. In this case, the operators A-, need 

not to be ordered according to their norms, so one 
can choose for A\ a non-hermitian operator. For ex¬ 
ample, in the case of the one-qubit decaying channel, 


^ dec — ^ 


= J 0 Vp 
0 0 


j, where p 6 [ 0 , 1 ] is 


1 0 

0 V^P 

a free parameter, the dual map <F[[ ec = {A\,Al} is 
not stochastic. However the inverted map & R ec = 

0 0 " r — 


VP 0 


V^P 0 
0 1 


]}• 


is stochastic. While the de¬ 


caying channel , F c jec describes the process of a sponta¬ 
neous decay ’downwards’ |1) —> |0) with probability p, 
the reversed process § R ec describes the transition ’up¬ 
wards’ | 0 ) —> | 1 ). Observe that for such a definition of 
the reversed map $ R ec the Crooks relation holds as a tau¬ 
tology. Note also that the invariants states of the map <F 
and its reverse are different. This is not the case for 
the time reversal operation $> Rc from [26| . Furthermore, 
in the case p > 0 the invariant state of the map <Fdec is 
pure and thus not invertible, so the operation is not 
well defined. 

Reversing Quantum Brownian motion We now revert 
back to the environmental representation (U) and assume 
that the system actually is connected to second physical 
system which acts as a heat bath. Such a map can be 
given by ((4|) where the ancilla state a is a thermal equi¬ 
librium state of a set of harmonic oscillators at inverse 
temperature /3, and U is a unitary time development of 
the system plus the bath determined by a total Hamilto¬ 
nian H = Hs + Hi + Hb where Hs is the system Hamil¬ 
tonian, Hb is the bath Hamiltonian and Hi is a linear 
interaction of the system and the bath. Standard time 
reversal of such a map is then given by (0, a procedure 
that can here be described as “attach the time-inversed 
bath, and run time backwards”. 

The quantum Brownian motion model is based on the 
observation that if bath frequencies form a continuum 
with Ohmic spectrum and a spectral cut-off, and the 
temperature goes to infinity, then © represents classical 
Kramers-Langevin dynamics p = /( x, t ) — r]jj + ^/2rj/j3f 
and x = -jj [H, H3] • Applying ([5j to the quantum Brow¬ 
nian motion model all operators are time reversed ac¬ 
cording to their parity, which in the classical limit means 
t t* = tf — t, Xt —> x'fY = Xt * and pt —> p/ + = 
— Pt *, and the Kramers-Langevin equation transforms 

int0 If = /0*> */ - **) + Vw + \/ 2 V/Pt an d %r = Ti- 

Time reversals of classical stochastic differential equa- 
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tions have been extensively discussed in the literature, 
and it is well understood that they are not unique E2- 
The example just derived by taking the classical limit of 
quantum Brownian motion is the “natural time reversal” 
of Kramers-Langevin dynamics, but also other possibili¬ 
ties make sense (note the “anti-friction”!). 

A second example of time reversal of Kramers- 
Langevin dynamics, in E2 called “canonical time rever¬ 
sal” , is based on the same variable transformation but 
assumes that the conservative force and the friction force 
transform differently under time inversion, resulting in 
a Kramers-Langevin dynamics also for the time-reversed 
motion i.e. aL p * = - t*) - Vjj + V 2 r )/P€ and 

= jj. To lift this definition to "R" we obviously 
have to treat the system and the bath differently. For 
the bath time must run forwards, so as to result in dis¬ 
sipation, while the conservative effects embodied in the 
force / are to be applied in the opposite order. This 
can be achieved by considering the system Hamiltonian 
of the form Hs = Hkin + V(xs,t) and the two unitary 
operators 

jj _ -y- e _ x(/o / H kin +V(x s ,t)+Hr+HB^ (16) 

JjRc ‘~TC~ ^ Hkin + V(. x S,tf— + + (17) 

where R stands for time ordering. Time reversal by 
changing U to U R amounts to the procedure of “attach 
the bath, let time run forwards, but time-reverse the ex¬ 
ternal drive”. In [l7j several other examples are given of 
time reversals of stochastic dynamics which can also be 


“lifted to R”. 

Discussion: In this work we have introduced a general 
notion of time reversals of quantum maps which general¬ 
izes standard time inversion in quantum mechanics. As 
in classical dynamics in contact with a heat bath, this 
definition of time reversal is not unique. One possibility 
- but not the only possibility - is to choose an envi¬ 
ronmental representation of the quantum map, and then 
apply standard quantum time inversion on the combined 
system and ancilla. Another possibility is to start from 
an intrinsic definition of the quantum map in terms of 
Kraus operators, and then define time reversal on that 
level. In any case, from any such definition one can de¬ 
fine an entropy production in the environment functional 
analogously to the classical setting, and for each such def¬ 
inition quantum fluctuation relations are satisfied identi¬ 
cally. 
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